Matter Wave Turbulence: Beyond Kinetic Scaling 
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Turbulent scaling phenomena are studied in an ultracold Bose gas away from thermal equilibrium. 
Fixed points of the dynamical evolution are characterized in terms of universal scaling exponents 
of correlation functions. The scaling behavior is determined analytically in the framework of quan- 
tum field theory, using a nonperturbative approximation of the two-particle irreducible effective 
action. While perturbative Kolmogorov scaling is recovered at higher energies, scaling solutions 
with anomalously large exponents arise in the infrared regime of the turbulence spectrum. The 
extraordinary enhancement in the momentum dependence of long-range correlations could be ex- 
perimentally accessible in dilute ultracold atomic gases. Such experiments have the potential to 
provide insight into dynamical phenomena directly relevant also in other present-day focus areas 
like heavy-ion collisions and early-universe cosmology. 
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I. INTRODUCTION 

Turbulence is a characteristic phenomenon naturally 
observed in many dynamical settings far from equilib- 
rium. Besides its generic appearance in the context of 
fluid dynamics, wave turbulence has become an impor- 
tant topic in non-linear dynamics. In laboratory exper- 
iments of ultracold atoms, recently turbulent dynamics 
of vortical motion in Bosc-Einstcin condensates has re- 
ceived increasing attention from both theory and experi- 
ment, see Refs. P, [HH, HI- Similar experiments may pro- 
vide insight into dynamical phenomena relevant also in 
other present-day focus areas like collision experiments 
of heavy nuclei and reheating after inflation in early- 
universe cosmology, see Refs. [a, S 0] for recent reviews. 

The theoretical framework for turbulence was formed, 
to a large part, in the 1940s Q. In the usual approach, 
wave turbulence is characterized by scaling laws derived 
in the framework of kinetic theory, see Refs. @, Q for 
reviews with further references. There are two require- 
ments for the validity of kinetic theory: The de Broglic 
wavelengths of particles must be small compared to the 
mean free path between collisions. This allows one to 
describe them as classical colliding particles. The second 
requirement is that the duration of individual scattering 
events should be small compared to the mean free time 
between collisions. Otherwise, interference between suc- 
cessive collisions would spoil their mutual independence. 

We study, in the present work, turbulent dynamics of 
an ultracold Bose gas beyond quantum kinetic theory. 
We use quantum field theoretical methods, with the aim 
of forming a theory framework for the systematic study 
of turbulence phenomena in experiments with ultracold 
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atoms. Non-thermal scaling solutions for the dynamics 
are found, and these are characterized, in the infrared, by 
anomalously large exponents, confirming results recently 
presented in Refs. [l(J [H| m t nc context of relativis- 
tic field theory. Such an extraordinary enhancement in 
the low-energy distribution of matter-wave modes could 
be experimentally accessible in dilute ultracold atomic 
gases. 

We point out in this work the strong connections be- 
tween the phenomenon of " weak wave turbulence" , which 
is perturbative and well described within kinetic theory, 
and the predicted non-perturbative low-energy scaling 
solutions. Since the latter occur in the deep infrared, the 
scaling behavior comes with distinctive properties such as 
critical slowing down and universality characteristic for 
critical phenomena. We explain in detail why this far- 
from-equilibrium critical dynamics may be understood 
as a " strong turbulence" phenomenon. As in the kinetic 
theory of weak wave turbulence, two different solutions 
are found in both the ultraviolet and infrared limits. It 
is shown that the obtained solutions corresponding to a 
quasiparticlc cascade in the ultraviolet regime require a 
fixed dispersion relation between frequency and momen- 
tum while it is sufficient to have energy conservation for 
the solutions which relate to a scale- invariant energy flux, 
i.e., energy cascade. 

Turbulence can arise as a stationary phenomenon in 
a driven system away from thermal equilibrium, as well 
as during a transient period after the driving force has 
ceased. It then offers the system a special path of transit 
to eventual thcrmalization. Progress in understanding 
the structure of state space with respect to trajectories 
connecting non-equilibrium with equilibrium has poten- 
tial applications in very different areas of physics. The 
preparation of ultracold atomic Bose and Fermi gases in 
various trapping environments allows to precisely study 
quantum many-body dynamics of strongly correlated sys- 
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terns. A number of experiments have recently focused on 
far-from-equilibrium dynamics, long-time evolution, and 
thermalization in such systems, cf., e.g., Refs. [H, EH, E3 • 
A possible realization of turbulent dynamics in an 
ultracold atomic gas could be achieved with a Bosc- 
Einstcin condensate mode as a source of excitations, 
which starts oscillating parametrically after a quench 
of the scattering length and thus the chemical poten- 
tial. The associated phenomenon of parametric reso- 
nance leads to a strong enhancement of infrared modes. 
This has been studied in great detail in weakly coupled 
theories of preheatin g a fter inflation in the early uni- 
verse EE EH E3, EM Ell- In this case, weak turbu- 
lence describes above a characteristic momentum scale 
the gradual transfer of energy towards higher momenta 
in the form of a cascade. 

In contrast, for low momenta occupation numbers 
grow non-perturbatively large and exhibit a different 
scaling behavior. A correct description of these low- 
energy modes requires non-perturbative approximations. 
A non-perturbative approach that has gained substan- 
tial success in describing far-from-equilibrium evolution 
over long times is based on the 2PI effective action 
in next-to-leading order (NLO)l/7V approxi- 
, |24| . Dynamic equations derived from an nPI 
effective action by definition obey crucial conservation 
laws, as for the total energy and, in a non-relativistic 
system, for particle number. Here we use the 2PI ap- 
proach in NLO 1 /Af approximation to study turbulence 
in ultracold gases. 

This article is organised as follows: Sect. [TT] provides 
the set of non-equilibrium dynamic equations which are 
used to find a necessary condition for the existence of 
a non-thermal stationary point. In Sect. IIII( we derive 
the stationarity condition and discuss the thermal fixed 
point and summarize results for turbulent scaling found 
within kinetic theory as far as relevant for the later dis- 
cussion. Sect. IIVI provides the derivation of the scaling 
behavior, both in the perturbative regime as well in the 
non-perturbative infrared limit of low-energy mode exci- 
tations. Our conclusions are drawn in Sect. Ivl Technical 
details and a short summary of the real-time functional 
quantum field theory are provided in the appendix. 




scattering between quasiparticlcs with a well-defined dis- 
persion relation. After successful applications of these 
nonperturbative expansions to the study of far-from- 
equilibrium dynamics and thermalization in relativistic 
bosonic [H, [il, [H, Hi and fermionic [13, H, Ell theo- 
ries, they have recently been employed in the context of 
ultracold bosonic quantum gases For 



introductory texts see, e.g., Refs. [35|, [36[. Some details 
about the 2PI effective action approach to nonequilib- 
rium dynamics relevant for our discussion can also be 
found in App. [Bl 

We consider the evolution of an ultracold bosonic 
quantum many body system described by the com- 
plex A/"-component Heisenberg field operators <I> Q (i,x), 
a = 1, . . . , Af in. d spatial dimensions, obeying the com- 
mutation relations [$ Q (t, x), y)] = 5 a p5(x — y), 
[$ a (t, x), <&p(t, y)] = 0. In the following we choose a 
basis where the field is written in terms of its real and 
imaginary components, $ Q = ($i, Q +i$2,a)/v / 2. Includ- 
ing the field component index as well as the "magnetic" 
index a into a single index a = (i a , a), with i a = 1, 2, the 
commutation relations read 



[$ a (t,x.),$ b (t,y)] = -a\ 



<5 a(3 (5(x 



y), (i) 



where a 1 denotes the Pauli 2-matrix [4^]. We consider 
a quantum field theory for a complex A/"-component field 
tp a (x) (a = (i a ,a), i a = 1,2, a = 1,...,JV) with quartic 
interactions, 



S[<p] 



1 



AAf 



(p a (x)iDj(x,y)(p b (y) 



fa(x)iPa(x)(f b (x)ip b {x), (2) 



where we use the notation J = J dxo J d d x with 
(xq,x) = (i,x). This model describes, e.g., an ultra- 
cold Bose gas of atoms with Af hyperfine sublevels whose 
interaction strength g docs not depend on the particu- 
lar hyperfine scattering channel of a pair of atoms. In 
d = 3 dimensions the coupling strength in such a system 
is g = Aira/m, a being the s-wave scattering length. The 
free inverse classical propagator of this model reads 



II. NONEQUILIBRIUM DYNAMICS BEYOND 
KINETIC THEORY 



In this section we summarize the aspects of nonequilib- 
rium quantum field theory most relevant to the discussion 
of turbulence in this article. Our analysis is based on the 
two-particle irreducible (2PI) effective action [2(| Hi], [22j 
applied to real-time dynamics. We approximate the 2PI 
effective action to next-to-leading (NLO) in the nonper- 
turbative expansion in powers of the inverse number Af 
of internal degrees of freedom [23|, [24| . This description 
can reach substantially beyond a quantum kinetic ap- 
proach, which is typically restricted to on-cncrgy-shcll 



<t>=0 



5(x- y)S a/ 3 [-i&i a i b d Xo - H 1B (x)S ab ] 



(3) 



see Eq. (|B5|) . Here H 1B (x) 



HU^/2m + V{x) 



denotes the single-particle Hamiltonian, and we choose, 
in the following, the external potential V(x) to vanish. 



A. Dynamic equations 

In this section we recall the equations which describe 
the time evolution of the lowest order connected corrcla- 
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tion functions or cumulants [33l . l35j | 

0«(ar) = (*a(i)>, (4) 
G ab (x, y) = (T$ o (a;)*&(l0) - M^Mv)- (5) 

The two-time correlation function G defined in Eq. (JSJ) 
involves a time-ordered product of in general noncom- 
muting field operators. Therefore, G can be decomposed 
as [13 



Gab(x, y) = F ab (x, y) - -sgn(xo - yo)p ab {x, y), (6) 

where the signum function sgn(xo — yo) evaluates to 1 
( — 1) for xq later (earlier) than yo, and where the statis- 



tical component F and the spectral part p are defined 
in terms of the anticommutator and commutator of the 
fields, respectively, 



F ab (x,y) = ±({$ a (x),$ b (y)}) c 
p ab (x,y) = i([$a{x),$ b (y)}}. 



(7) 
(8) 



Here, (-) c is a short-hand notation for the cumulant 
($ a $ b ) c = ($ Q $ h ) - ($„}<$6>, see Eq. ©. 

The resulting integro-differential dynamic equations 
for 4>, F, and p read for Gaussian initial conditions 



( - i(j2 abdxo ~ gF ab (x,x)^j(j} b (x) - (h 1b (x) + ^[(f> c (x)(/) c (x) + F cc (x,x)]^(t)a(x) = ( dy^ ib (x,y;<j) = 0)(j) b (y), (9) 



[^oo^xo +M ac (x)] F cb (x,y) 
[i°ic d xo + M ac{x)\ Pcb(x, y) - 



J.'o 

to 
■I'd 



dzE* c (x,z;(/>)F cb (z,y)+ / dz Ef c (x, z; 4>)p cb {z, y), 



dz^ c (x, z; 4>)p cb (z,y). 



(10) 



(11) 



These identities are equivalent to Kadanoff-Baym or Schwingcr-Dyson equations. Here we employ the notation 
J. dz = J. dzo J d d z. The "mass" matrix M containing the free Hamiltonian as well as mean-field potential terms 



is defined as 



M ab (x) = S ab H 1B {x) + !j-(<f> c (x)(f> c (x) + F cc (x,x)j +gU>a (s) 4>b (x) + F ab (x,x) 

I 



(12) 



The self energy E acts as the kernel in the non-Markovian 
memory integrals in the above dynamic equations and ac- 
counts for collisions building up correlations in the sys- 
tem. It is obtained as the derivative of the 2PI part I^, 

^ab\x, y; (/>, G) — 2i— — -, (13) 

5G a b{x,y) 

and has been decomposed into a local mean-field part 
T/i){x) adding to the mass matrix, and a nonlocal part 
written in terms of statistical and spectral components, 

^ab(x,y) = E%}(x)6(x - y) 

+ ^ab{x,y) ~ ^sga(x -y )H p ab {x,y). (14) 

These non-local parts form the kernels for the memory 
integrals on the right-hand sides of the integro-differential 
dynamic equations (f9|)- (fTTj) . 

B. NLO 2PI 1/JV expansion 

To practically solve the dynamic equations (f9"l)- (fTTj) . 
details about the self energy S arc required, and these 



are, in general, only available to a certain approxima- 
tion. In the following we employ an expansion of I^ in 
powers of the inverse number of field degrees of freedom 
TV [HIM OH. The l/M expansion to next-to- leading or- 
der (NLO) is equivalent to replacing certain vertices in a 
loop expansion by a bubble-resummed vertex [M HH • In 
the context of an ultracold Bose gas, it has been discussed 
in Refs. [13, HI, HI, H3| ■ This approximation scheme has 
also been recovered in a functional rcnormalization group 
inspired approach [38j where it results as a truncation 
in orders of proper n-point functions combined with an 
s-channel approximation of the equation for the proper 
four-vertex. 

In this scheme, the contribution T2\<f>, G] to the 2PI ef- 
fective action involves a leading (LO) and next-to-leading 
order (NLO) part which can be diagrammatically repre- 
sented as shown in Fig. [1] in terms of 2PI closed loop 
diagrams involving only bare vertices, full propagators 
G, and field insertions <j>. While the leading-order contri- 
bution involves one diagram, in NLO a chain of bubble 
diagrams is resummed. All of these diagrams are propor- 
tional to the same power of 1 /AT since each vertex scales 
with 1/Af, which is cancelled by the (blue) propagator 
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r 2 LO [<j),G] 



r 2 NLO [<j>,G] : 



o o 



Or * /v-O-t* 



O 



-x 



O 



FIG. 1: (Color online) Diagrammatic representation of the 
leading order (LO) and next-to-leading order (NLO) contri- 
butions in the l/A/'-expansion, to the 2PI part T2[4>, G] of the 
2PI effective action. The thick blue lines represent 2-point 
functions G a b{x,y), the red crosses field insertions <j> a (x), and 
the dashed lines vertices g/4J\fS(x — y). At each vertex, it is 
summed over double field indices a and integrated over double 
time and space variables x. 



loops which scale with Af since they involve a summa- 
tion over the field indices from 1 to Af [10, [H[ . 

Inserting the 2PI effective action expanded to NLO in 
l/Af into Eq. (fT3"| one calculates the self energy E and 
from this the local and nonlocal contributions defined 
in Eq. (|14[) . The resulting expressions are provided in 
Appendix 



III. STATIONARY SOLUTIONS OF THE 
EVOLUTION EQUATIONS 

In this section, we set up the condition for station- 
ary solutions of the dynamic equations introduced above 
which is later used to derive turbulent scaling behavior. 
We furthermore discuss the thermal equilibrium solutions 
as well as the scaling solutions associated to weak wave 
turbulence found in the framework of kinetic theory. 



A. Stationarity condition 

We consider stationary homogeneous solutions of the 
dynamic equations ([9])- (jlip . i.e., solutions invariant un- 
der translations in time and space, 

<t>a (x) = <Pa(t) ~ cxp(-i/it) ; 
F a b(x,y) = F ab (x - y); p a b{x,y) = p a b{x - y) (15) 

Note that the stationarity condition for cf> allows for a 
remaining rotating phase with constant angular veloc- 
ity corresponding to a chemical potential, i.e., it only 
requires the condensate density \<p\ 2 to be constant in 
time. 

As we are looking for stationary solutions, we can send 
the initial time in the dynamic equations to to — > — oo. 

F o 

Consequently, also the self-energy components £ a £ , as 
well as the functions I F,P and P F ' P which are defined in 
Eqs. (|C4p - (IC8[) depend only on the relative coordinate, 



^ab( x 'V) = ^ab( x ~ v)> etc - We wr ite the Fourier trans- 
form of the statistical correlation function as F a i,(p) := 
J d d+1 x exp(ipx) F a b{x) , etc. for all other two-point func- 
tions. 

A necessary condition for a solution of the dynamic 
equation (|10p for the statistical correlation function F 
to be translationally invariant as defined in Eq. (I15|) 
reads 



J(p) 



Kb(P) F ba(p) 



^ab(p)Pba{p) = 0. 



(16) 



In Appendix [DT] we derive the necessity of condition |T| 
for stationarity of F, using the symmetry properties of 
the correlation functions involved. This stationarity con- 
dition can also be obtained from the trace of a gradient 
expansion of the evolution equation (|10|) . which leads to 
2u)dx F(X ,Lj, p) = J(Xo,u,p), with u = po [H]. In 
situations where F and p are independent of the time 
Xq = {xq + yo)/2, the gradient expansion in Xq involved 
in the derivation of this equation becomes exact. We re- 
mark that there is no analogous condition following from 
the dynamic equation pip for p. Following the same ar- 
guments as for F, we find that Eq. (fTTj) allows for trans- 
lationally invariant solutions obeying the conditions (|15p , 
see Appendix ID II 



B. Thermal equilibrium 

Before proceeding to nonthermal stationary solutions 
we verify that in thermal equilibrium the condition (|16p 
is satisfied: In this limit, the grand canonical density 
matrix reads p = exp(—/3(H — pN)/trexp[—f3(H—pN)j, 
and one can deduce the fluctuation-dissipation relation, 



F ab\u,p) = -i ( UBe(u} 



Pcb ( W >P) 



with 



"-BE 



M = l/(e^-") - 1), 



(17) 



(18) 



uj = p Q . For a detailed discussion see, e.g., [3a, [36|: In 
addition to this relation there is a similar relation for the 
self-energies, 



(w,p) = -i [ jibeM + ^ 



(19) 



Substituting this equation and Eq. (jTTJ) into Eq. (fl6|) one 
finds, as expected, in thermal equilibrium the stationarity 
condition is fulfilled. 



C. Weak wave turbulence from kinetic equations 

In this subsection we briefly review the Kolmogorov 
theory of turbulent scaling solutions of the wave kinetic 
equation Q. Going away from thermal equilibrium, a 
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kinetic description of the time evolution implies a quasi- 
particle assumption, which assumes a fixed dispersion 
relation ui = w(p) between momentum and frequency. 
Similar to (JT7J) one writes for spatially homogeneous sit- 
uations, with t = X°, 

F <« p >( W ,p) = -i (n^(u) + 1) Pa f\u;,p), (20) 

where the quantity ri/ qp ) (uj) will play the role of a quasi- 
particle number of mode p for fixed uj = ut(p). The 
spectral function is given, e.g., for a cold Bose gas in the 
symmetric, i.e., non-condensed phase, with quadratic dis- 
persion m(p) = p 2 /2m, by the ideal-gas expression (cou- 
pling g = 0) [3J] 

* ) (p) = pS ) (p)=-(%-^) + ^ + ^)) 

x [sgn(w) S ab ~ o- 2 ab ], (21) 

where a 2 is the Pauli 2-matrix. 

Inserting this ansatz into the dynamic equation (|10p 
for F, where the self energy has been approximated by 
the second-order-in-g expressions given in Appendix [C] 
one derives, in leading order of a gradient expansion in 
t = Xq 7 and for <p = 0, the quantum four- wave kinetic 
equation [34l [36} 

d t n p = I(p,t), (22) 

I(P, t) = g 2 J d d k d d q |T pkqr | 2 <5(p + k - q - r) 

X S(0J p + Ld k - CJ q - W r ) 

x [(n p + l)(n k + l)n q n r 
- n p n k (n q + l)(n r + 1)], (23) 

where n p = n^(p,t). In the case defined by Eq. (f2~Tjl . 
the transition matrix clement squared |Tp kqr | 2 is a nu- 
merical constant independent of momenta. The pertur- 
bative expansion for weak coupling g underlying (|23p 
restricts the occupation numbers to be parametrically 
n p <C (|p|a) _1 . For 1 <C n p -C (|p|a) _1 , where a is 
the s-wave scattering length, the equation goes over to 
the classical kinetic equation. This is the regime, where 
scaling solutions describing Kolmogorov wave turbulence 
can be obtained. In the quantum limit n p <C 1, waves 
behave like particles, and Eq. (]22[) reduces to the Boltz- 
mann equation. This is typically the range of high mo- 
menta, where occupation numbers are low, and no tur- 
bulent scaling will be observed in the quantum regime. 
Owing to the local conservation of n p , the above kinetic 
equation can be written as a continuity equation 

dtn^ p \p,t)+d t Q 1 (p 1 t) =0, (24) 

with the divergence of the current, diQi(p, t) = —I(p, t), 
defined in terms of the scattering integral /, Eq. (f2"3"l) . 

Alternatively one can write the kinetic equation (1221) 
as a continuity equation for the energy density e(p, t) = 



u(p)n( qp ' (p, t), with the divergence of the current P of 
energy flow given by the integral diPi(p, t) = — I(p, t) = 
-w(p)J(p,t). 

Taking into account the spherical symmetry of the in- 
teractions, the scattering integral I(p, t) can be aver- 
aged over the spatial directions which returns an integral 
I(uj(p)) depending on the frequency only. As is shown in 
detail in Ref. 0, the stationarity condition I(uj(p)) = 
has four universal scaling solutions (w(p) = o>(|p|)) for 
the classical kinetic equation: 

n (qp) (w(sp)) = n( qp )(s|p|) = S - K n( qp )(|p|), (25) 

where the exponent is either k = and k = z, for the 
constant and thermal solutions, respectively, z being the 
scaling exponent of the quasiparticle frequency uj(sp) = 
s z uj(p), or 

k = kq = —(3d + 2m — z), (26) 
3 

■y 

K = Kp = KQ + -. (27) 

Here, m is the scaling exponent of the transition matrix 
element, T(sp, sk, sq, sr) = ,s m T(p, k, q, r) = s m Tp kqr , 
which evaluates to m = for the weakly interacting 
cold Bosc gas away from unitarity. The above nonther- 
mal scaling exponents characterize the Kolmogorov sta- 
tionary distributions of wave turbulence. The solution 
with exponent kq corresponds to a |p|-indcpendent ra- 
dial particle flux |p| d_1 Q(|p|), while the solution with 
exponent Kp gives a momentum- independent energy flux 
|p| d_1 P(|p|) Q. For the non-condensed cold Bose gas, 
one finds the momentum scaling exponents 

KO = d-|, (28) 
k p = d. (29) 

For a system with linear dispersion, z = 1 and dominant 
four-wave interaction with m = — 2z = —2, e.g., an ul- 
tracold Bose gas in the unitary limit at large momenta 
|p| (neglecting inelastic scatteri ng) , or a relativistic scalar 
theory in the high-energy limit [361 ] , one recovers the scal- 
ing exponents 

5 

K Q = d--, (30) 
4 

K P = d--, (31) 
which in d = 3 evaluate to kq = 4/3 and kp = 5/3 [46j]. 

D. Stationarity condition for the full dynamic 
equation 

We proceed by considering the stationarity condition 
(|16p beyond the above kinetic approximation, i.e., using 
the l/Af expansion of the 2PI effective action to NLO 
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+ 




2PI 1/N: 

X - X ■ > < + >0 < 

FIG. 2: (Color online) Next-to-leading order (NLO) contri- 
butions to the scattering integral J. Upper panel: Diagram- 
matic representation of the NLO contributions in the 1/M- 
expansion, to the self energy £ a &(a;, y; <j>). The big filled circle 
represents a resummed vertex which in NLO 1/Af is defined 
as shown in the lower panel. The (red) (j> dependence of the 
vertex in the two-loop diagram contributing to E indicates 
the further internal field dependence integrated over within 
the loops as can be seen from applying Eq. (|13p to the lower 
line of diagrams of r£ LO in Fig. [1] All other symbols are as 
in Fig. Q] 



following Rcfs. A major restriction of the ki- 

netic description is that it cannot describe high occupa- 
tion numbers n p beyond the range of validity of pertur- 
bation theory. Since scaling solutions n p ~ |p|~ K imply 
large occupancies at low momenta, kinetic theory breaks 
down in the infrared. In contrast, the non-perturbativc 
2PI 1/Af expansion to NLO allows us to investigate the 
low-momentum regime, which is characterized by differ- 
ent scaling solutions than the above Kolmogorov results. 

The NLO 1/Af self energy entering J(p), Eq. (gO)), is 
depicted in Fig. [2j Before proceeding to the scaling be- 
havior of J we derive, in this section, the formal expres- 
sion for J in terms of the correlation functions F and p. 



1. Vanishing field expectation value 

We first consider the simpler case of a vanishing field 
expectation value <p a = 0. Scaling will result for strongly 
occupied modes where quantum fluctuations can be ne- 
glected. Expressed in terms of the statistical and spectral 
components of G this means that F a b(x) 2 ^> p a b(x) 2 . 
In thermal equilibrium this is readily obvious from the 
fluctuation-dissipation relation (|17p . As was shown in 
Rcfs. [H, H3| away from equilibrium the classical statis- 
tical limit implies that the p ■ p terms are neglected in 
the function I F defined in Eq. (|C5j) and, for = 0, the 
I p ■ p term in E F , Eq. (|C4[) . The translationally invariant 
self-energies then read 



where 



~X(I F -F ab )(x) 

-\[(I P -F ab )(x) + (I F -pab)(x)], (32) 



A = 2g/Af, 



and 

I F = (1 - I R ) * U F + I F * IT 4 , 

p = (i-i R )*n p + i p *u A , 

with the retarded and advanced functions 



(34) 



I R {x) = (0 ■ I p )(x), I A {x) = (6- ■ I")(x), (35) 
U R (x) = (6 ■ U p )(x), Il A (x) = {0- ■ Il p )(x). (36) 

In the classical limit the functions H F,P arc 

n F = A(F ■ F)/2, U p ^X(F-p). (37) 

Here and in the following we use the more convenient 
notation in terms of products and convolutions of func- 
tions in x as well as p space, see App. [Alfor details. The 
product/convolution of two correlation functions implies 
sums over field indices, (F-p)(x) = F a (,(x) -p a b{x). Finite 
integration limits in time in the above convolutions are 
taken into account by the theta function 6(x) = 0(xq), 
with 9~ (x) = 6(—x). Overall arguments (x) have been 
suppressed. 

The recursive equations (|34|) for I F and I p can be 
solved explicitly. Using the convolution theorem one ob- 
tains, in momentum space, 



i f (p) = (\ cS -n F )( P ), 

I p ( P ) = (\ cS -W)(p), 



(38) 



where Il F (p) = \(F * F)(p)/2, W(p) = X(F * p)(p), and 
1-I R 1 



A eff (p) 



i-n- 4 (i -n A )(i + n R ) 
l 



\i + m 



(39) 



Here, we have suppressed arguments (p). The second 
equality in Eq. ([39| follows from 



(l+U R (p))-(l~I R (p))=l 
which in .T-space reads: 



(40) 



I1 R - I R - U R * I R = 0. (41) 



This identity is proven by substituting the expression 
([34]l for I p into the second term giving 



n J 



e-[n p -(9-i p )*w+i p ' 



■n p ))-n R *(0-i p ) = o 

(42) 

which is verified using (|36|) as well as that for any two 
functions / and g: ((0 ■ /) * (6 ■ g)) = 6 ■ ((0 ■ f) * (0 ■ g)). 

We are now ready to write down an expression for J 
in the case of vanishing <fi: Substituting Eqs. (|38p into 
(33) (|3"2"j) and the self energies into condition (fT5|) . we obtain 
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for J» = J^=°(p): 



reads 



J°(p) = -\{Fba ■ [F ab * (2A cff • (F * p)) 
+ Pat * (A cff • (F * F))] 
- Pba ■ [F ab * (A cff • (F * F))] } 

\2 f 

= ; — ^ d d+1 kd d+1 qd d+1 rS(p + k-q-r) 

2(2tt) 8 J y ' 

x A off (p + k) [F ab (p)F ab (k)F cd (q)p cd (r) 

+ F ab (p)F ab (k)p cd (q)F cd (r) 

- F ab {p)p ab {k)F cd {q)F cd {r) 

- p ab {p)F ab {k)F cd (q)F cd {r)]. (43) 



2. Non-vanishing field 

In the general case one needs to take into account the 
possibility of a non- vanishing condensate field <p 7^ 0. 
As pointed out above, scaling will require that statis- 
tical fluctuations dominate over quantum fluctuations, 
F ab (x) 2 ^> p ab (x) 2 . In this case the self-energies (IC4|) 
read 



Ef 6 (a:) 
KbW 



-,F(0=O) 

,p(0=O) 
J a6 



(44) 



A (I^afo + P P ■ F ab + P F ■ p ab ) 



(45) 



J X (P) = 



A 



J\ P ) 



A :! 



^a0 6 ■ [Pah ■ (F * F) ■ A cff 

+ 2 Pcd ■ (F dc * (F afc • A cff )) 

- 2F ab ■ (F * p) ■ A cff 

- 2F cd • (F dc * (p ab ■ A cff )) 

- 2F cd .(p dc *(F afc -A off ))], 

[pba ■ {Fab * [A ■ (F * F)}} 

2F ba -{F ab *[{F*p)-A]} 

Fba ■ {Pab * [(F * F) ■ A]}, 



with the effective coupling 

A(p) = (f>a(/> b A ab ■ A off . 



(49) 



(50) 



(51) 



J x and J A written in the form of momentum integrals 
as in Eq. (|43[) for J° are provided in App. ID 21 To make 
contact with the notation of Ref . [l(| , we introduce 



J 4 = J° + J A and J 3 = 



(52) 



The integral J 4 results from the left self-energy diagram 
in the upper panel of Fig. [21 the integral J 3 from the right 
one in which the field <f> is not integrated over inside the 
loops. This concludes the calculation of J. We have three 
terms J = J° + J x + J A given by Eqs. (|43|). (1D10|) . and 

flnnj) . 



IV. TURBULENT DYNAMICAL SCALING 



Here and in the following we suppress arguments, e.g. 
I F = I F (x), where they can be inferred. 

I F and I p are as before, and P F , P p are rewritten 
analogously as described in Appendix ID 2[ giving 

P F (p)=\cl> a <t> b [F a b-\ eS ~I F -A ab ] (46) 
P"(p) = \^ b [ Pab ■ A off - I p ■ A afc ] , (47) 



where (see Appendix ID 2 



A ab (p) = m 



1 + U R 



(48) 



We are now ready to derive the full function J. The 
ingredients I F , I p , P F , and P p are given by Eqs. ([38]) . 
(|4^)) . and (|47p. respectively. Substituting these into the 
full self-energies, and (|55)). we get J = J° + where 
J° was obtained in Eq. (|4"3")) and the field dependent part 



A. The scaling ansatz 

We are looking for solutions that fulfill the stationarity 
condition (fTB]) in the infrared. Following Refs. PUSH we 
consider scaling solutions with properties 

Pab(s z po, sp) = s~ 2+v p ab (p 0l p) (53) 
F ab {s z p a ,sp) = s- 2 - K F ab (p 0l p), s>0. (54) 



The free spectral function (|2lj) fulfills the scaling prop- 
erty (j53]) . with z = 2 and 77 = 0. Deviations from this 
scaling arc accounted for by a modified value for z and 
an "anomalous dimension" r\ (see e.g. Ref. [4l|) different 
from zero, which we will assume to be small in the follow- 
ing. The scaling ansatz (|54|) for F is chosen analogously, 
with an exponent k still to be determined. 

The scaling relation ([53)) . with 77 = 0, remains valid 
beyond the case of an ideal gas, for the near-zero- 
temperature weakly interacting Bose gas in the bro- 
ken phase which is described by Bogoliubov theory and 
merely a free gas of quasiparticlcs interacting with the 
condensate mode. In this case, the dispersion becomes 
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linear in |p| at low momentum, reflecting the sound- wave 
character of the modes and requiring z = 1 in the scaling 
relations, but quadratic scaling is ensured by the Bogoli- 
ubov coefficients. At zero temperature one has 

Pabip) = in [8(po - ^ P ) + $(po + w p )] 

x [(ul + vl)6 ab sgn(p ) -al b ], (55) 

where uj p = [e p (e p + 2g<j) a (j) a )} 1 / 2 (with e p = p 2 /2m) 
is the Bogoliubov dispersion which is linear in |p| at 
p 2 <C ig(j) a (j) a . u p = [(e p + g(j) a 4>a + uj p )/2uj p ] 1/2 and 
v p = I( e p + g^'ai'a — u p )/2tu p ] 1 ' 2 are the Bogoliubov co- 
efficients resulting from diagonalization of the propagator 
and scale like u p ~ |p| -1 ^ 2 , v p ~ | p> | 1 / 2 in the infrared. 
Hence, the scaling (|53[) also applies to function ([551) . 

In the limit of large occupation numbers, the density of 
particles with momentum p in momentum space, n(p) = 
(tft(p)#(p)>, is given by 

n(p) = \J^ \Fn(p) + iF M - iF 2 i(p) + F 22 (p)} . 

(56) 

Hence, the scaling behavior (jM)) of F implies that 

ra(sp) = s 2 ^ 2 " K 7i(p). (57) 

Since in addition we assume isotropy, it follows that 
n(p) = n(|p|) ~ |p| z ~ 2 ~ K . Hence, for a quadratic single- 
particle dispersion, z = 2, k acts as the occupation- 
number scaling exponent. 

In order to derive the scaling exponent k from the nec- 
essary condition (|16[) for stationarity we take into account 
the spatial isotropy of the Hamiltonian and consider the 
ensuing weaker necessary condition which results from 
integrating Eq. Q16p over the spatial momenta p [47| . 

J( Po ) ee J d'p [E p ab (p)F ba (p) - ^ b (p) Pba (p)] = 0. (58) 

In the following we determine, under the assump- 
tion that the correlation functions obey the above scal- 
ing properties, for which values of k the stationarity 
condition (f58|) is fulfilled. To derive these values we 
make use of scaling transformations [13, fill ] which im- 
ply reparametrizations of the frequency-momentum inte- 
grals in J(j>o), similar to the Zakharov transformations 
used in the context of weak wave turbulence 0] . The thus 
reparamctrized integrals can be transformed to a unique 
form with the help of the scaling laws for the correla- 
tion functions which relate these functions at different 
frequency and momentum scales to each other. 

In the following we will consider the scaling behavior 
in the ultraviolet (UV) and infrared (IR) regimes sepa- 
rately. We first derive the scaling exponents in the UV 
regime. They correspond the perturbative Kolmogorov 
exponents, k p , Eq. ([25]). and kq, Eq. (|25jl . We will then 
present the derivation of a different scaling regime, with 
significantly enlarged values for exponents, in the IR. 



B. Scaling exponents in the ultraviolet regime 

The aim of this section is to show that in the perturba- 
tive regime of sufficiently large frequencies and momenta 
the exponents characterizing weak wave turbulence 

8 n 

K = d+z-- + i (59) 

and 

« = d+|(*-2) + | (60) 

allow J(j)q), Eq. (|58p . in d spatial dimensions to van- 
ish. Here Eq. (|5_9|) corresponds to Eq. j28]) and Eq. ([60]) 
to Eq. (|2"9")l for z = 2 and 7] = 0. The conditions for 
the existence of these scaling solutions arc different. The 
solution with k given by Eq. (|59p requires that the ex- 
citations of the system are described sufficiently well by 
quasiparticles, Eq. (|20p. with a fixed dispersion relation 
(f2Tj) while the solution (|60|) exists also without this re- 
striction. In the UV the leading-order perturbative ap- 
proximation with a zero- width spectral function becomes 
valid such that kinetic theory and therefore both scaling 
solutions apply. 



1. Dominant contribution to J 

J has three components J = J° + J x + J A , given by 
Eqs. (O, IjDlOp . and (|D11|) . We derive scaling laws for 
J°, J A and J A in the UV limit to determine which of 
the contributions to J(j>) dominates in this regime. For 
sufficiently large frequencies and momenta the occupa- 
tion numbers are low. As a consequence, one can classify 
the contributions to Eq. (|58p perturbatively according to 
powers of A. Since any non-zero width of the spectral 
function is of higher order in A, a (5-like spectral func- 
tion can be assumed under the integral of Eq. (|55|) to a 
given order in the coupling. A similar reasoning applies 
to the perturbative statistical function. Therefore, we 
will recover the results of kinetic theory. 

Turning first to J° we need to find the scaling law 
for the effective coupling A eff defined in Eq. (|39|) . As 
discussed in Appendix [El scales like H R (s z po, sp) = 
s z+d - 4 - K+ m R (p ,p) (see Eq. fE4])). If k> z + d-A + rj 
which according to Eq. (|60p requires r\ < z/2 + 2, this 
implies that |n fl (p)| < 1 in the high- momentum limit. 
Hence, one can neglect in the denominator of (|39p . 
and obtains: 

A off (p) w 1 (UV limit). (61) 

Away from the UV limit the above relation becomes ap- 
proximative as does the scaling behavior discussed in 
the following. This also applies to the kinetic wave- 
turbulence theory where one needs to show the local- 
ity of the scattering integral. Locality implies, that the 
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scaling law at a particular momentum scale is indepen- 
dent of what happens far away in momentum space. In 
kinetic theory, the locality of the scattering integral fol- 
lows from the convergence of the integral Q. We note 
that a proof of convergence of the scattering integral in 
the full dynamical theory is more difficult than in the 
kinetic approximation, where the dispersion is fixed and 
the frequency integrals can be eliminated. The existence 
of approximate scaling solutions at finite momenta can 
be shown by means of numerical simulations of the dy- 
namics as presented, e.g., in Refs. [6LHoj]. 

It follows that in the UV limit, we can neglect the ef- 
fective coupling in J°(p), Eq. (|43|) . The same arguments 
apply to the remaining contributions J A and J A , and in- 
serting the scaling pr opert ies (|53jl, ([54]) . p6]) . (|E9]) . and 
(|ETT|) into Eqs. (gSJ) , jUTO)) . and (|PTT|) one finds the scal- 
ing relations 

J°(s z p , sp) = s^ d+z -^-^J°( P(h p), (62) 
JVpo, sp) = s d+z - 6+ "- 2K J x ( Pa ,p), (63) 
J Vpo, sp) = s^ d+z - 5 ^- 3K J A (p , p). (64) 

We anticipate the general result k — d + 4(z — 2)/3 + 
?7/3. Comparing the scaling laws for J°, J x and J A , 
which scale with exponents 2(d+z — 4) +77 — 3k = — d— 2z, 
d + z-6 + i]-2K = -d- (5z + 2 + ?/)/3 and 2(d + z-5) + 
2rj — 3k = —d— 2z — 2 + 77, we come to the conclusion that 
both J° and J A dominate in the limit of large momenta 
for z = 2, assuming 77 = 0, while for z < 2 or 77 > 
the integral J°, and for z > 2, the integral J A dominate 
on their own. It will turn out that the same value of 
« = d + 77/3 renders J° and J A to vanish simultaneously 
if z = 2. 



2. Scaling transformations 

Taking into account the possible dominant contribu- 
tions in the UV regime (see last paragraph) the station- 
arity condition (|58p for the momentum integral over J(p) 
requires 



J°(po) = 0. 



(65) 



We derive a relation between k, 77, and z which fulfills 
this condition. Using the above argument that we can 
set A eff (p + k) = 1 in the integrand, the integral reads 



J°(po) 



A 2 



I d d+1 kd d+1 qd d+1 * 



2(2tt) 8 

x 8{p + k — q — r) 

x [Fab{p)F ab (k)F cd (q)p cd (r) 
+ F ab (p)F ab (k)p cd (q)F cd (r) 
- F ab {p)p ab {k)F cd {q)F cd {r) 
~ Pab(p)F ab (k)F cd (q)F cd (r)]. 



As is described in more detail in App. [F] it is possible 
with the help of scaling transformations to exchange p 
with another integration variable in Eq. ((6"6"|) , thereby 
keeping po as a free variable [h], HH • These transforma- 
tions are similar to the Zakharov transformations typi- 
cally employed to compute exponents for weak wave tur- 
bulence for a fixed dispersion relation and z = 2, T) = 0. 
In that case, frequency integrals are separated into differ- 
ent domains of integration, which can be mapped onto 
each other Q. Here we use scaling transformations to 
map the different integrands of the integral Eq. (|66jl . 
Since we do not assume a fixed dispersion relation in this 
approach, we will be able to use the same method also 
to compute IR scaling exponents in the non-perturbative 
regime below. This transformation makes use of the scal- 
ing properties (|5"3"|) and (f5"4")) . As these scaling relations 
involve an s z in the po component, one first rewrites 
Eq. (|66[) such that the integration variables fco, qo , and 
r are positive: 



dfc d<7odr /(jj, k, q, r) 



dfc d<7odr 



fco>0,g >0,r >0 



x (f(p, k, q, r) + [r -> -r] + [q -> -q] + [k -> -k] 



[q- 


> -q, r - 


> — r] + [k — > —k, r — » 


—r 


[k- 


*-k,q- 


> -q] 




[k- 


*-k,q- 


>-q>r-> -r]). 





(67) 

Here, the function / contains the full integrand in 
Eq. (|66p . including the remaining spatial integrals. For 
Po > 0, the fifth summand in paranthescs vanishes due 
to the delta- function S(p + k — q — r), such that there are 
4-7 = 28 terms contributing to the integrand in Eq. (|66p . 
Negative arguments can easily be dispelled using the 
symmetry properties for F and p: F ab (~p) = F ba (p), 

Pab{-P) = -Pba(p)- 

At this point we use a scaling transformation (see 
App. [F| to achieve that in each of the 28 summands, p 
carries the p-argument. We furthermore use that we can 
permute the integration variables k, q, r at will. Com- 
bining terms that are equivalent one arrives at the final 
result: 



d^ / d d+1 /fcd d+ V d+1 7 



(66) 



J°(po) = 

&o>0,g >0,r >0 

ab ab cd cd a a a a a \ 

X (S p+ k- q -r PpFkFqFr Po (Po + k Q ~ W ~ r o) 

ab ba cd cd —A A A A A 

+ 5 p -k+q+r PpFkFqFr Po (Po ~ K +% + r ) 

ab ba cd cd _ » AAA A\ 
+ 5 p - k - q - r PpFkFqFr Po {.Po - k ~% - r Q ) 

ab ab cd dc _a, A A A A\ 

+ 28pJ r k-q+r PpFkFqFr Po \Po + ^0 

ab ba cd dc a a a a a \ 

+ 25 p - k - q+r PpFkFqFr Po (Po - fco - 9o + r o )). 

(68) 

ab 

Here we have used the short-hand notation p p = p a b{p), 
cf. Eq. (|D12[) . and 5 p = 5(p). Combining the scaling of 
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correlation functions, delta function, and integral mea- 
sures, the exponent A results as 

A = -(3k - 3d + 8 - rf) - 3. (69) 

We now consider the condition that J x (po), Eq. (|D10|) 
vanishes, 

J X (po) = 0. (70) 

As before, if po is chosen large, only the high-momentum 
behavior of the functions involved is important, in par- 
ticular A cff ?s 1. Rewriting the integral such that all fre- 
quency variables are positive and neglecting terms van- 
ishing by the delta function, one obtains 9 • 3 = 27 sum- 
mands. By use of a scaling transformation of the form de- 
fined in App. [Fl Eq. (|F3|) . each of the 27 terms is rewrit- 
ten such that p carries the p argument. Analogous steps 
as before lead to the final expression 

J X (p )=<l> a <t>b / d d p [ d d+1 kd d+1 q 





fc >0,<3o>0 








X (Sp-k-q 


ab cd cd 
PpFkFq 


-A/ A 

Po (Po 


- k A 



-?o A ) 


+ Sp-k+q 


ab cd dc 
PpFkFq 


-A/ A 

■Po (Po 


Z-A 




+ 4) 


+ Sp+k-q 


ab cd dc 
PpFkFq 


-A/ A 

■Po (Po 


+ fco A 


-4) 


+ ^Sp-k-q 


cd dc ab 
P pFkFq 


-A/ A 

■Po (Po 


- k A 

K o 


-4) 


+ 26 p -k+q 


cd dc ab 
PpFkFq 


-A/ A 

■Po (Po 


- k A 




+ 'Zfip+k-q 


cd cd ab 
PpFkFq 


-A, A 

■Po (Po 


+ k£ 


-1o)), (71) 



with 

A= -(2«;-2d+6-7?) -2. (72) 



3. Scaling exponent: Results in the UV regime 

The delta distributions ensure that J°(po), Eq. (|68[). 
vanishes if 

A = l <^ K = d+ -(z-2) + ^. (73) 
J x (po), Eq. (TTTj) . becomes manifestly zero if 

A = l ^ K = d+|(z-2) + |. (74) 

In the special case of a vanishing field expectation value, 
ip = 0, J = J , and the result for k remains unaf- 
fected. We therefore come to the conclusion that in a 
cold Bose gas with quadratic dispersion, z = 2. and i] = 
Kolmogorov-likc scaling occurs in the UV regime, with 



k = d. The exponent ([60)) implies a scaling of the parti- 
cle number, Eq. (|57|) . as 

n(|p|) - |p|- d -(-- 2 +'))/3. (75) 

The scaling exponent k = d obtained for z — 2 and ?/ = 
corresponds to the Kolmogorov exponent Kp, Eq. I|29[). 
for weak wave turbulence in the kinetic approximation, 
see Sect. IIII CI We now consider the derivation of kq, 
Eq. (EH). 

Inserting the perturbative quasiparticle behavior, 
Eq. ((20|) , into Eq. (|68|) causes the first four lines in round 
parentheses to vanish because, after evaluating the fre- 
quency integrals over fco, qo, and ro, and thus setting 
the frequencies to the respective quasiparticle frequen- 
cies, qo = w(q), etc., one finds that 

^(wk.kjFaj^k'.k') = f ab (o;k,k)F ab (a; k ^k') = 0, 

(76) 

for arbitrary momenta k, k', while in general 
Pab (u k ,-k)F ba (uv,V) + 0, ^(wk.k^^.k') + 0. 
As a consequence, only the last term cx p A — k A — q A + r A 
remains in Eq. (|68|) . and J x (ui p ) vanishes identically, 
cf. Eq. (|7Tj) . Hence, a further scaling exponent results 
since J°(j>o) vanishes now also for 

A = & K = d + z-| + J. (77) 

For z = 2 and r\ = we obtain n = d — 2/3. Wc empha- 
size that a fixed dispersion relation or (|20p is required to 
find scaling with this exponent. Without (j20"|) . A = 1 
leaves in general nonvanishing contributions to J (po), 
Eq. Ij55]). as well as to J x (p ), Eq. 1(71]). These terms are 
in general nonzero since the frequency dependences of F 
and p allow for collision events which, as is read off the 
delta functions in, e.g., Eq. (|68| . change the number of 
wave excitations. Hence, the solution Kg, Eq. (f2"9")l only 
holds in the quasiparticle limit, Eq. (|20|) . This finding 
is consistent with the fact mentioned in Sect. IIII CI that 
Kolmogorov's exponent kq results under the assumption 
of a p- independent radial quasiparticle flux Q(|p|) which 
in turn requires the existence of a well-defined quasipar- 
ticle number. 



C. Scaling exponents in the infrared regime 

In this section we show that in the non-perturbativc 
regime of small frequencies and momenta a different scal- 
ing behavior arises than what is found above for weak 
wave turbulence. The exponents 

K = d + z-T) (78) 

or 

K = d + 2z-r/ (79) 
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allow J(po), Eq. (|58p to vanish in the infrared. Despite 
the fact that there are strong corrections to kinetic the- 
ory in this regime, we point out that certain properties 
of the non-perturbative solutions Eq. ([75]) and Eq. ([75]) 
are still similar to the Kolmogorov solutions Eq. (|59|) 
and Eq. (|60p . respectively. Major differences concern, 
apart from the numerical values, properties such as criti- 
cal slowing down or universality for the infrared behavior, 
which is characteristic for critical phenomena, here far 
from equilibrium. Universality of a critical phenomenon 
requires that the integral equations, which determine the 
critical exponents Eq. (|75|) and Eq. (f?9"]) . are dominated 
by the low-frequency and -momentum limit such that 
higher momenta do not affect the value of exponents. 
We also emphasize that, while the perturbative UV so- 
lutions are only valid at high momenta if quantum cor- 
rections are neglected, the non-perturbative IR solutions 
are a property of the quantum theory. Stated differently, 
classical-statistical fluctuations dominate over quantum 
corrections at low momenta, where occupations numbers 
are high, such that the quantum and the classical theory 
are characterized by exactly the same exponents. 

1. Dominant contribution to J 

J has three components J = J° + J A + J A , given by 
Eqs. ([33"]). (|D10p . and (|Dllj) . Again, we first derive scal- 
ing laws for J°, J x and J A in the IR limit to deter- 
mine which of the contributions to J(j>) dominates in 
this regime. Turning first to J° we need to find the scal- 
ing behavior of the effective coupling A cff . As discussed 
in Appendix [El if k > d + z — 4 + 77, which according 
to Eq. ((7JJ) requires r\ < z/2 + 2, then > 1 in 

the denominator of A off , Eq. (|39|) . in the low- momentum 
limit. Hence, A cff , in the IR limit, scales like 

X cS ( s z Po , S p) = s 2(K+4 - 2 - d -'">A cff (p ,p) 

= s 2 ( z+4 - 2 ")A eff (p ,p). (80) 

In the last line we anticipated the result ([79]). Now 
consider J°(p) in the low-momentum range. Proceed- 
ing through the same steps as for the UV case one finds 
the IR scaling of the three integrals contributing to J(p) : 

J°( S z Po, S p)= S - K -' 1 J a (p ,p), (81) 
J x (s*p , sp) = S ~ d - z+2 -^J x ( Po , p), (82) 
J A ( S ^o, sp) = s- d - z+2 -iJ A (p Q , p). (83) 

Comparing these scaling relations we conclude that for 
k > d + z — 2, the integral J° dominates in the IR limit. 
Given the result ([79]) this condition requires rj < z + 2. 

2. Scaling transformations 

Taking into account the dominant contributions in the 
IR regime (see last paragraph) the stationarity condition 



([58| reads 

•A;Po)=0. (84) 

In the following we show that this leads to a different 
value for k as compared to the UV regime. We obtain, 
for J°(po), an expression identical to Eq. ([66]) . except 
that an overall factor A (p + k) multiplies the integrand, 
cf. Eq. (|43p. It is a remarkable property that the main 
consequence in the non-perturbative regime is the ap- 
pearance of a momentum-dependent effective coupling. 
As a consequence, in this sense there is a direct link be- 
tween the weak-wave-turbulence analysis in the UV and 
the calculation of exponents in the infrared regime of 
what may be called strong turbulence. 

We proceed like in the UV case, splitting the integral 
such that it contains only positive frequencies, applying 
scaling transformations of the type (|F3[) and finally com- 
bining terms. We obtain a result similar to {68]), but 
a different expression for A: As A now contributes a 
scaling factor this results as (compare Eq. ([69]) ): 

A = -{k + t] - d- z) (85) 

z 

Hence, J°(po) vanishes identically if 

A = l n = d + 2z-T]. (86) 

This result also holds for the special case of vanishing 
field, (f> = 0, where J = J°. The exponent ([79]) implies a 
scaling of the particle number, Eq. ([57]) . as 

n(|p|) ~ | P r d -*- 2+ ". (87) 

A second solution is obtained from Eq. ([85]) for A = 0. 

A = n = d+z- 1]. (88) 

Given our discussion in Sect. IIV B 31 this solution re- 
quires a fixed dispersion relation. While in the UV this 
is described in terms of a perturbative, i.e. J-like spectral 
function, this is not the case for the infrared properties of 
critical phenomena. Here the power-law behavior of the 
spectral function leads to dominant contributions in the 
infrared, i.e. for both p° and |p| approaching zero. For a 
detailed discussion of this in the context of equilibrium 
critical phenomena, see Ref. [42]. In Refs. [|| [l5| numer- 
ical simulations of the classical equations of motion for 
a relativistic scalar theory with quartic self-interactions 
were presented which demonstrated the evolution of the 
system into a turbulent scaling regime after parametric 
resonance, confirming perturbative results for the expo- 
nents of weak wave turbulence. New numerical simula- 
tions were presented in Ref. [Io| . which extend to the 
infrared and demonstrate the presence of a strong turbu- 
lence regime with strongly enhanced correlations. Using 
the 2PI effective action in NLO 1/Af approximation, the 
properties of this scaling behavior in the infrared regime 
are recovered analytically (Io| . The simulations indicate 
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that after parametric resonance or spinodal decomposi- 
tion dynamics, the infrared scaling exponent k = d+z—r] 
is approached taking z = 1 and r\ = 0. We emphasize 
that this is the exponent which requires a fixed disper- 
sion relation as it is not expected in the IR regime. Hence 
the numerical results of Ref. [Icj . together with our re- 
sults indicate that the finite width of the dominant peak 
of the spectral function in the IR limit p — > allows for a 
scaling solution reminiscent of turbulence, corresponding 
to the Kolmogorov scaling with kq in the UV. More- 
over, according to the numerical results this should be 
approximately valid and representing a form of "strong 
turbulence" at small finite momenta. 



D. Thermal scaling 

We close our analysis by considering the scaling at the 
thermal fixed point. The obtained exponents k = d+z—r) 
and k = d+z — 8/'3 + ri/3 in the IR and UV limits, respec- 
tively, are larger than the exponent for thermal equilib- 
rium: From the fluctuation-dissipation relation (|17p and 
([T5)l . it follows that, in the Raleigh- Jeans limit jibe 3> 1, 
i.e. f3u ~ and hence e 13 ^ »1 + j3u), that 

F ab (u, p) ~ (Puj)- 1 ■ p ab (u, p) (89) 

and hence 

F ab (s z u,sp) = S - 2 - z+r >F ab {Lo,v). (90) 

As a consequence, the thermal scaling exponent in the 
large-momentum limit reads 

k = z — t] (thermal equilibrium). (91) 

Note that in a Bose gas in the broken, i.e., condensate 
phase, at low momenta k is modified by the interac- 
tions. Using the Bogoliubov dispersion cj p = [e P (e p + 
2g4> a (f) a )] 1 / 2 , with e p = p 2 /2m, which becomes linear in 
|p| at p 2 <C 4gcf) a (j) a one finds that, at zero temperature 
and with r\ = 0, 

k=1 (IR limit, T = 0). (92) 

Taking into account Eq. (|57[) . n = z — r/ leaves, however, 
the scaling of the particle number invariant as compared 
to the UV, scaling, 

»(|p|) ~ |Pr 2+ "- (93) 

We note that it was shown in Ref. (43j . using sum rules 
for moments derived from linear response theory that 
2n p + 1 > 2kBTm<f> a (f> a /p' 2 , which implies that, at finite 
temperature, 2 — z + k > 2, i.e., 

k>z (IR limit, T > 0). (94) 

in the infrared limit. 



V. CONCLUSIONS 

We have presented a scaling analysis of wave-turbulent 
fixed points in the dynamical evolution of a spatially uni- 
form ultracold Bose gas away from thermal equilibrium. 
The analysis was focused on the scattering (self-energy) 
term of the dynamical equations for two-point correla- 
tion functions, in nonperturbative next-to-leading-order 
approximation of a l/Af approximation. The dynamic 
equations in this approximation were derived from the 
2PI effective action which ensure the approximated time 
evolution to obey crucial conservation laws as those for 
total energy and particle number. 

Searching for a dynamical fixed point was performed 
by implementing the necessary condition that the mo- 
mentum integral over the scattering term vanishes in a 
nonthermal stationary state. Presupposing scaling prop- 
erties for the correlation functions entering this condi- 
tion, scaling exponents were derived for the two-point 
functions characterizing a nonthermal stationary state. 
Such a state, showing nontrivial scaling properties, is 
well known in the theory of wave turbulence. We repro- 
duce the perturbative Kolmogorov scaling properties dis- 
cussed, e.g., in Ref. 0, identifying it as the scaling behav- 
ior in the high-frequency regime. We find that the tur- 
bulent spectrum, i.e., the momentum mode occupation, 
in d spatial dimensions, behaves like \p\^ d ~( z ~ 2+r v/ 3 m 
this regime, where d is the number of spatial dimensions, 
\p\ z the scaling of the dispersion, and rj an anomalous 
dimension assumed to be small. In deriving this scaling, 
a possible nonzero constant condensate field was taken 
into account, but the two-to-two scattering events inde- 
pendent of this field were found to dominate the scaling. 

As in the kinetic theory of weak wave turbulence, two 
different solutions are found in both the ultraviolet and 
infrared limits. Solutions which correspond to a quasi- 
particlc cascade in the ultraviolet limit require a fixed dis- 
persion relation between frequency and momentum while 
it is sufficient to have energy conservation for the set of 
solutions which relate to a scale-invariant energy flux, 
i.e., an energy cascade. In the UV the leading-order per- 
turbative approximation with a zero- width spectral func- 
tion becomes valid such that kinetic theory and there- 
fore both scaling solutions apply. At momentum scales 
approaching zero the divergence in the spectral function 
dominates its behavior such that also in the IR limit both 
scaling solutions are relevant. 

At low momenta and frequencies we find strong correc- 
tions to the kinetic theory underlying the perturbative 
Kolmogorov scaling analysis. A different scaling regime 
appears in the infrared, with significantly enhanced scal- 
ing exponents resulting from the nonperturbative char- 
acter of the interactions in the gas. Here the nonper- 
turbative nature of the 2PI l/JV expansion to NLO is 
essential to be able to describe this physics. Our findings 
confirm analogous results presented for the relativistic 
case in Refs. [[(J El- This scaling enhancement should 
show up in a strong occupation of low-momentum modes, 
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rising with a power law n(p) ~ |p|~ d ~ 2 ~ 2+? ', as com- 
pared to the above quoted Kolmogorov scaling for high 
momentum modes. Since the phenomenon occurs in the 
deep infrared, the scaling behavior comes with distinctive 
properties such as critical slowing down and universality 
characteristic for critical phenomena. Our analysis shows 
that this far-from-equilibrium critical dynamics can be 
understood as a strong turbulence phenomenon. 

Based on our results we propose to study turbulence 
phenomena in dilute ultracold gases with hindsight to the 
scaling properties of correlation functions, in particular 
of the momentum distribution of particles in the gas. It 
would be striking to find experimentally the predicted 
strong turbulence regime, which would have important 
consequences also in other areas of physics such as early- 
universe cosmology where similar phenomena can cru- 
cially determine the thermal history of our universe. 
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The convolution theorem is then 
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APPENDIX A: NOTATION 

We chose the (H ) convention for the Minkowski 

metric. The Minkowski product between 4-vectors 
p = (p°,Jj\25 2 ,J5 3 ) = (po,p) = (w,p) and x = 
(x°, x , x 2 , x 3 ) = (xo,x) = (£, x) is then px = poXo — 
p • x. Functions of 4-vectors are Fourier-transformed 
with the Minkowski product as follows: T[f(p)](x) = 
f(x) = (2n)-*fd 4 pexp{-ipx}f(p), ^[/(x)](p)_ = 
f(p) = J d 4 xexp{ipx}/(x). In d < 3 spatial dimensions 
all definitions arc analogous. 

Function letters denote both the function and its 
Fourier transform. The argument specifies which of the 
two is meant. The following convention is used for con- 
volutions: 



(f*g)(x)= / d d+1 yf(y)g(x-y), 



(f*g)(p) = 



d d+1 q 



f(i)g(p- <?)• 



(Al) 
(A2) 



APPENDIX B: THE 2PI EFFECTIVE ACTION 
APPROACH TO NONEQUILIBRIUM 
DYNAMICS 

In this appendix we review the aspects of the 2PI ef- 
fective action approach to nonequilibrium dynamics. For 
introductory texts see, e.g., Refs. [H, [36J. 



1. 2PI effective action 

For a given initial-state density matrix po{to) charac- 
terizing a system far from equilibrium, all information 
about the quantum field theory is contained in the gen- 
crating functional for correlation functions: 

Z[J, K; p D ] = Tr [p D (t ) T c cxp {»( J ^ (x)$ a (x) 

y a (x)R c ab (x, y )<f b (y))}], (Bi) 



4>„ 



xy,C 

with Heisenberg field operators < & a (x) obeying the non- 
relativistic commutation relations (JTJ) . In Eq. (|B1[) . Tq 
denotes time-ordering along the closed time path C lead- 
ing from the initial time to along the real time axis 
to some arbitrary time t and back to to, with J x c = 

J c dxo J d d x. Contour time ordering along this path cor- 
responds to usual time ordering along the forward piece 
C + and antitemporal ordering on the backward piece C~ . 
Note that any time on C~ is considered later than any 
time on C + . The source terms in Eq. (|Blj) allow to gen- 
erate correlation functions by functional differentiation 
such as 



(T c $(x 1 )---$(x„)> = 



5 n Z[J,K; PD ] 



i n 5J{xi) ■ ■■8J{x n ) 



J.K=0 



(B2) 
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where the field indices have been suppressed, Xi = 
(a;o,i,Xi) is a d + 1-dimcnsional coordinate, and we have 
used that for the closed time path Z = 1 in the absence 
of sources. 

In order to obtain the dynamical evolution of the sys- 
tem, i.e., the time dependence of the correlation functions 
(|B2[) requires precise knowledge of the generating func- 
tional Z . For most nontrivial practical applications it is, 
however, not feasible to directly compute Z , e.g., by use 
of Monte Carlo techniques. In contrast to many cases of 
imaginary-time evolution in equilibrium this is prevented 
by a variant of the sign problem. A possible way out is 
the reformulation of the problem in terms of an effec- 
tive action functional. The final aim of this approach 
is the derivation of time evolution equations for certain 
correlation functions of interest, in particular for the low- 
est order connected correlation functions or cumulants as 
defined in Eqs. (gj) and ©. 

The functional derivative relations (|B2[) which deter- 
mine these two correlation functions are rewritten into 
stationarity requirements of an effective action functional 
r[<^>, G] with respect to its arguments <f> and G. V is ob- 
tained by a double Legendre transform of the Schwingcr 
functional W[J, R] = —i\nZ[J,R] with respect to the 
source fields J and R introduced in the definition (jB2|) 
of the generating functional Z . The result is the 2PI 
effective action 



r[^G] = r( llo °p)[0,G] + r 2 [0,G], 



(B3) 



where 



r (lloop)^ G ] = + ^ h n Q- 



G 1 G) 



const. 



(B4) 

denotes the one-loop part obtained in a saddle-point ap- 
proximation of the generating functional Z involving the 
classical action S and the classical inverse propagator 



5 2 S[4>] 



8(t>a{x)8(t>b{y) 



(B5) 



The trace in Eq. (|B4[) includes sums over all internal 
and space-time indices. The functional S depends on the 
classical field only and defines the underlying model to 
be considered. A particular class of models considered in 
this article will be defined in Sect. Ill Al below. 

The stationarity requirements for T with respect to <fi 
and G, 

T~TTT =~Ja(x)- [ R a b{x,y)(j) h {y), (B6) 



smc] i 

77T-7 r = -^Rab{x,y), 

dG ab {x,y) 2 



(B7) 



are equivalent to the equations of motion of 4> a and G a b- 
If the initial state po{to) is taken to be Gaussian, 
then only the above introduced cumulants a (i o ,x) and 
Gab(to-, x; to, y) are nonzero at t = to. Then, the ini- 
tial state pz)(to), in the generating functional Z, can be 



absorbed into the sources J and R at t — to. As a con- 
sequence, these sources can be set identically to zero in 
the equations of motion (|B6[) , (|B7p which thus close and 
describe the evolution of an isolated system starting in a 
state defined only by the first and second order cumulants 
at initial time to- 



2. Dynamic equations 

The effective-action approach defined in Eqs. (|B3[) . 
(|B4|) is, so far, valid for a general class of models de- 
scribing the evolution of a nonrclativistic bosonic field 
$ a (x). In Eqs. we have specified the particular 

model by defining the classical action functional S. 

Given the classical action, the contribution r 2 [(/>,G] to 
the 2PI effective action, see Eq. (|B3[) . is defined in terms 
of an infinite series of all possible 2PI diagrams without 
external legs, formed from the bare 4-vertex defined by 
the interaction part in Eq. @, full propagators G a b{x, y), 
and classical field insertions cf) a (x), i.e., of all such dia- 
grams which do not fall apart on opening at most two 
lines G. From the resulting 2PI effective action the dy- 
namic equations are determined by functional differenti- 
ation as defined in Eqs. (|B6p and (|B7[) . Since in general 
r 2 is an infinite series, also the resulting equations con- 
tain an infinite number of terms of increasing order in the 
number of bare couplings, internal lines G and fields <p. 
To obtain a practically solvable set of equations the dia- 
grammatic series r 2 needs to be truncated as discussed 
in more detail in Sect. HlBl in the main text. 

While a more detailed account of the standard proce- 
dures summarized above is beyond the scope of this ar- 
ticle we point out that the principal reason for deriving 
the equations of motion for <fi and G by use of a station- 
arity condition is that the resulting equations conserve 
particle number and energy irrespective of the particular 
truncation of r 2 . In any such truncation the equations of 
motion imply a totally conserved particle current density 
j K (x) = (p(x), j(x)) as well as energy momentum tensor 
T K \(x) with d K j K (x) = and d K T K \(x) = 0, respectively. 
Integrated over space, these relations describe the con- 
servation of total particle number and energy. This im- 
portant feature forms a necessary requirement for the 
resulting equations to be applicable to long-time dynam- 
ics, i.e., most notably, not to lead to secular evolution in 
this limit. 



APPENDIX C: 2PI SELF ENERGY 

From the 2PI effective action the self energies (fl4|) in 
two-loop order, for cj> = 0, are found as 



hKb( x >y) 



29 



R ab( X ' v) 

hKbi^y) 
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H F (x,y) \W{x,y) \ ( F ab (x,y) 
■\TP(x,y) n F (x,y) J {-±p ab (x,y) 



P"(x,y;I F <n = -^{H p (x,y) 



with (F 2 = F ab F ab , etc.) 



TL F ( Xl y)\_ g ( F{x,yf - \p{x,yf 



and 



n p (z,y) J Af I 2F ab (x,y) Pab (x,y) 



R p 



2g(F 3 - \(Fp 2 + pFp + p 2 F) 



Af \ F z p + FpF + pF 2 



1 „3 



(CI) 



(C2) 



(C3) 



where (F 3 ) ab (x,y) = F ac (x,y)F dc (x,y)F db (x,y), (Fp 2 ) ab 
(x,y) = F ac (x,y)p dc (x,y)F db (x,y), etc. 

For the 2PI effective action in NLO of the expansion 
in 1/Af, T 2 [4>, G] = T%°[(f>, G] + T^ LO [(j), G] the diagram- 
matic expansion of which is shown in Fig. [1] we obtain, 
using Eq. (fT3|) , the self energies (|T4")) as 



Z F b (x,y)\ = _2 l 



I F (x,y) 



A F (x,y) i\AP(x,y) 
■\M(x,y) A F (x,y) 



\l p (x,y) ' <^(aO<M2/) 

F a b(x,y) 
\pab{x,y) 



(C4) 



where A F 'P(x,y) = I F >P{x,y) + P f >p(x, y; I F 'P). The re- 
summation to NLO in 1 /Af is taken care of by the coupled 
integral equations for I F ' P [3(| : 



I F (x,y)\ fU F (x,y) 
IP(x,y) J \U p (x,y) 

to dzlP ^ Z \w{ Zl y) 



U F (z,y) 



The functions P F ' P , which contribute to A F ' P in the self 
energies (|C4[) and vanish if 0j = 0, read [33[ 



AT 

dz [H F (x, z)IP(z, y) + I F (x, z)HP(z, y)] 
dz [HP(x, z)I F (z, y) + P(x, z)H F (z, y)\ 



to 

xo rVo 

dv / dwI p {x,v)H F (v,w)I p (w,y) 

to J to 

dv j dwI p (x,v)H p (v,w)I F (w,y) 

to J to 

fVo ryo 

I dv 

to J v 



I 1 
dwI F {x,v)H p (v,w)I p (w,y)[, (C6) 



x 

dz [HP(x, z)IP(z, y) + I p (x, z)H p (z, y)\ 

Vo 

%o fVo , 

dv / dwI p {x,v)H p {v 1 w)I p {w,y)\, (C7) 

' Vo J Vo 

wherein the functions H F,P arc defined as 

H F {x,y) = -<f> a (x)F ab (x,y)<f> b (y), 

H p (x, y) = -<j> a {x) Pab (x, y)4>b{y). (C8) 



APPENDIX D: STATIONARITY CONDITION 

1. General discussion 

In this appendix we present a proof of Eq. (|16[) : Con- 
sider the dynamic equation (|10p for translationally func- 
tions F, p, and 4>, obeying (fTS"|) . with y set to zero: 



[icr 2 ac d X0 + M ac ]F cb (x) 



d d+1 zEP ac (x-z)F cb (z) 



with 



+ / d d+ W ac {x-z)p cb {z) 



M ac = S ac [- J* + 9 -{4> d <t>d + F dd (0)) 
2m 2 

+ g(<Pa<t>c+F ac (0)). 



(Dl) 



(D2) 



Contracting over indices a = b, and using the symme- 
try properties a 2 ah = -of , F ab (x) = F ba (-x), F ab (x) = 
-F ba {-x) and (d 2 F ab )(x) = (d?F ba )(-x), one finds that 
the left-hand side of Eq. (|Dlj) is invariant under x — > —x. 
Considering the right-hand side of Eq. (|D1|) . and us- 
ing the further symmetry properties p ab {— x) = —p ba {x), 
V F ab {-x) = £f» and = -T, p a {x) following 

from Eqs. (|C4 |) -([C8 ]l one finds, adding Eq. (|DTjl and mi- 
nus its counterpart with the sign of x reversed and sum- 
ming over a — &, that 

= / d d+1 z {Z p ab (x - z)F ba (z) - T, F b (x - z) Pba {z)) 



(D3) 

which, by the convolution theorem, proves the stationar- 
ity condition (fl6|) . Following the above line of argument 
for the dynamic equation (|11[) for p immediately shows 
that the condition corresponding to Eq. (|D3|) is automat- 
ically fulfilled. 



2. Field dependent contribution to J 

In this appendix wc first sketch the derivation of 
Eqs. (|46[) and (|4T|) and finally quote the momentum- 
integral expressions for J A (|4^|) and J A (|50p . In the 
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translationally invariant case, cf. Eq. (|15[) . the integrals 
P F and P p can be written, using the notation introduced 
in App. \M and Sect. IIIIDli as 



P F (x) = A0 a b [F ab + F ab * I A + I F * G 



ab 



I" * Fab 

I A * F nb * I R 



+ I A *G A *I F ], 

\<t>a<Pb [Pab - Gf h * I" + p ab * I A 



: Pab 



G R b *I p *I R 

Pab *I A *I R 

G A *I A *I?], 



with the retarded and advanced propagators 



G R (x) = 9-p abl G A b (x) = 9- ■ Pab . 



(D4) 



(D5) 



(D6) 



The above expressions for P F and P p can be cast into a 
compact form with the help to A cff and A ab , to be defined 
below. Consider the complex conjugate of Eq. P0|) . 



(1 



P) ■ (1 - IP 4 ) = 1 



(D7) 



integral expressions for J A , Eq. (|49p . and J A , Eq. ([50)) : 

\2 

J*(p) = 



J A ( P ) 



A eff 

- a; 
K 



p 

eff 
fe 

eff 



ab cd cd ab cd cd ab cd cd 
PpFkFq ~ FpPkFq ~ FpFkPq 



cd ab dc 
PpFkFq 

cd dc ab 
PpFkFq 



cd ab dc 
FpPkFq 

cd dc ab 
FpPkFq 



cd ab dc 
FpFkPq 
cd dc ab 
FpFkPq 



(D10) 



2(2tt) 6 



d d+1 kd d+1 qd d+1 rd v+k - a - r A 



-\-k — q — r^p-\-k 



ab ab cd cd ab ab cd cd 
PpFkFqFr + FpPkFqFr 



ab ab cd cd 
FpFk P qF r 



ab ab cd cd 
FpFkFq P r 



(Dll) 



Here we have introduced a compact notation for matrix 
indices and momentum dependence, e.g. 



ab 

Pp= Pab{p)- 



(D12) 



Expression (|D11[) can be combined with J°, Eq. (|43|) , to 
J 4 , Eq. (52|) to 



2(2tt) 8 
x Xf +k • (1 - XMb^ab) 



' ab ab cd cd ab ab cd cd 
PpFkFqFr FpPkFqFr 



ab ab cd cd 
FpFk P qFr 



ab ab cd cd 
FpFkFq P r 



(D13) 



obtained with IL R (p)* = -Ii A {p), IP{p)* = -I p (p) and 
0(p)* = 0-{p). Hence 



\ eS (p) 



(i + n*)(i-n^) 



{i - e * i p ){\ + e- *ip). 

(D8) 



We furthermore define 



A ab (p) = 25? 



G 



1 + U R 



r<R 

U ab 



~ U ba 



i + u R i-n A 

* Pab) - (0~ * Pba) - {0 * Pab) ' (0 * I P ) 

- * Pba ) ■ (0~ * IP). (D9) 



Now, going over to momentum space and rearranging 
terms in Eqs. ()D4[) and (|D5[) . one obtains Eqs. (|46|) and 
63). 

We close this appendix by quoting the momcntum- 



APPENDIX E: SCALING PROPERTIES 

The scaling properties for p and F are given in 
Eqs. (|53|) and |54|) . respectively. In general, if two func- 
tions / and g scale like 

f(s z p 0l sp) = s~ 7 f(p) and g(s z p , sp) = s~ s g(p), 



it follows that their convolution scales like 



if * g)(s z p ,sp) 



z+d-~f-S 



(f*g)(p)- 



Hence, H p (p) = X(F * p), scales as 

UP(s z p ,sp) = s z+d - i+2 ' n - K UP(p), 

The Fourier-transformed step function 0{jp) 
J d d+1 x e ipx 9(x) scales as 

9( S z Po ,sp) = s- z - d 9(p). 



(El) 
(E2) 

(E3) 



By the scaling rule for a convolution (|E1[) . we see that 
the retarded function II (p) — 9 *H P scales like H p : 



U R {s z p ,sp) = s z+d - i+2 ^- K U R {p). 



(E4) 
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If k > z + d - 4 + 277, \U R (p) I > 1 in the IR, and one can 
neglect the 1 in the denominator of A cff = |1 + n^|~ 2 , 

A Gff (s^ ,sp) = s 2(K+4 -^' i - 2,?) A cff (p) (IR). (E5) 
In the UV limit, |II fl (p)| < 1, thus 

A off (» = 1 (UV limit). (E6) 



According to Eq. (|E1|) . the retarded propagator ()DGl) . 
G ab(p) = 0* Pab, scales like p: 

G R b (s z Po , S p) = s -^G R b ( P ), (E7) 

A ab (p) = 2R[G R b /(l + 11*)], Eq. jMl), as 

A ab ( s z Po , S p) = s K - d+2 - z -iA ab (p) (IR), (E8) 
A ab ( s z p ,sp) = S - 2+,? A ab (p) (UV). (E9) 

Finally, A(p) — (fafbAab ■ A c ^, Eq. (|51l) . scales as 

A( S 2 p ,sp) = s 3(«-d- 2 )+5(2-,) A ^ ^ E1Q ) 

A(a z po,sp) = S " 2+r 'A(p) (UV). (Ell) 

APPENDIX F: SCALING TRANSFORMATIONS 

Scaling transformations are extensively used within the 
Kolmogorov-Zakharov theory of wave turbulence Q . In 
this appendix we briefly describe the scaling transfor- 
mation used in Sects. IIVBI and IIV CI to derive scaling 
exponents. 

The scaling transformation allows to exchange two 
variables in an integral expression despite the fact that 
one of them is a free variable, if the functions of the in- 
tegrand obey a scaling property: Consider an integral 
which has the form of the scattering integral J 3 (p) inte- 
grated over the spatial momenta, J 3 (po) = J X (po), scc 



Eqs. (IDIOI) . IP]) (v= (p ,p), etc.): 

J(po)= J d d pd d+1 kd d+1 q5(p-k-q)f(p)g(k)h(q), 

k ,q >0 

(Fl) 

with po > 0, and /, g and h obeying the scaling laws 

f(s z p ,sp) = s" 7 /(p ,p), 
g(s z p ,sp) = s~ 5 g(p ,p), 
h(s z p ,sp) = s~ e h(p ,p). 



(F2) 



The following transformations allows to swap p and k [10l . 
El: 



po = j7 fc o; 

K 



, PO PO , 

ko = tjPo; qo = -rrio 



{ k) k; k ^UJ p; q= U' 



(F3) 



leading to the result 



J(po)= / d d pd d+1 kd d+1 q8(k-p-q)f(k) g {p)h(q) 

k o ,qo>0 

x ( Po /k o r P (F4) 

where 

- (3=-(2d + 2z-j-6-e). (F5) 

z 

Hence, the arguments p and k are exchanged at the cost 
of a factor (po/fco) ^ which involves the scaling expo- 
nents. The scaling transformations applied in Sects. HVBl 
and IIV Cl on the integrals contributing to J 4 (po) which 
involve one more frequency-momentum integration are 
performed analogously. 
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